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Abstract 

I present a more refined analysis of the mass-splittings between the different heavy quarkonia 
states, using new double ratios of exponential moments of different two-point functions. Then, I 
test the validity of the 1/m-expansion, extract a s (Mz) — 0.127±0.011 from M Xc ^ P \^— M^ps), and 
provide a new estimate of the gluon condensate from — M Vc and M Xb — Mx, which combined 
with the recent estimate from the r-like decay sum rules in e + e~ — > J = 1 hadrons data, leads 
to the update average value (a s G 2 ) = (7.1 ± 0.9) x 10~ 2 GeV 4 from the light and heavy quark 
systems. I also find Mr — M Vb ~ 63^5? MeV implying the possible observation of the rjb in the 
T-radiative decay. 



Talk given at QCD 9 6- Montpellier, 4- 12th July 1996 and at the 28th ICHEP96-Varsaw, 25-Slth 
July 1996, and based on the paper \hep-ph/95 1 234$ (to appear in Phys. Lett.B386, October 1996), 
which will be referred as SN. 
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1 The double ratio of moments 



QCD spectral sum rule (QSSR) a la SVZ M (for a recent review, see e.g. 0) has shown since 
15 years, its impressive ability for describing the complex phenomena of hadronic physics with 
the few universal "fundamental" parameters of the QCD Lagrangian (QCD coupling a s , quark 
masses and vacuum condensates built from the quarks and/or gluon fields). In the example of the 
two-point correlator: 

Tl Q (q 2 ) = tjd'x (0\TJ Q (x) (J Q (o)) t |0>, (1) 

associated to the generic hadronic current: Jq(x) = QTQ(x) of the heavy Q-quark (r is a Dirac 
matrix which specifies the hadron quantum numbers), the SVZ-expansion reads: 

D=0,2,... dimO=D (Mq — q 2 



(2) 

where « is an arbitrary scale that separates the long- and short- distance dynamics; are the 
Wilson coefficients calculable in perturbative QCD by means of Feynman diagrams techniques; 
(O) are the non-perturbative condensates of dimension D built from the quarks or/and gluon fields 
(D — corresponds to the case of the naive perturbative contribution). Owing to gauge invariance, 
the lowest dimension condensates that can be formed are the D = 4 light quark m q (ipijj) and gluon 
(a s G 2 ) ones, where the former is fixed by the pion PCAC relation, whilst the latter is known to 
be (0.07 ± 0.01) GeV 4 from more recent analysis of the light [[| quark systems. The validity 
of the SVZ-expansion can be understood, using renormalon techniques (absorption of the IR 
renormalon ambiguity into the definitions of the condensates and absence of some extra l/g 2 -terms 
not included in the OPE) j|, || and/or by building renormalization-invariant combinations of the 
condensates (appendix of || and references therein). The SVZ expansion is phenomenologically 
confirmed from (among other applications) the unexpected accurate determination of the QCD 
coupling a s @-|| and from a measurement of the condensates || from semi-inclusive r-decays. 
The previous QCD information is transmitted to the data through the spectral function IhiIIq^) 
via the Kallen-Lehmann dispersion relation (global duality) obeyed by the hadronic correlators, 



which can be improved from the uses of different versions of the sum rules |I],[0],H-|T2]. In 
this paper, we shall use the simple duality ansatz parametrization: "one narrow resonance" + 
ll QCD continuum" , from a threshold t c , which gives a good description of the spectral integral in 
the sum rule analysis, as has been tested successfully in the light-quark channel from the e + e~ — > 
1=1 hadron data and in the heavy-quark ones from the e + e _ —* i(j or T data. We shall work 
with the relativistic version of the Laplace or exponential sum rules: |1|, [H| [11], [12] : 



f°° 1 
Cni&jm 2 ) = / dt exp(— ta) — ImlWt), 

Jim? n 



Kh(o-) = --^\ogC H (o-,m 2 ), (3) 

where the QCD expression known to order a s is given (without expanding in 1/m) in terms of the 
pole mass m(p 2 = m 2 ), [] and contains the gluon condensate (a s G 2 ) correction to leading order. 
To this order, the gluon condensate is well-defined as the ambiguity only comes from higher order 
terms in a s , which have, however, a smaller numerical effect than the one from the error of the 
phenomenological estimate of the condensate, a = r = 1/M 2 (notations used in the literature) is 
the exponential Laplace sum rule variable; H specifies the hadronic channel studied. In principle, 
the pair (a, t c ) is free external parameters in the analysis, so that the optimal result should be 



^or consistency, we shall work with the too- loop order a s expression of the pole mass [O. 



insensitive to their variations. Stability criteria, which are equivalent to the variational method, 
state that the best results should be obtained at the minimas or at the inflexion points in a, while 
stability in t c is useful to control the sensitivity of the result in the changes of t c -values. These 
stability criteria are satisfied in the heavy quark channels studied here, as the continuum effect 



is negligible and does not exceed 1% of the ground state contribution [Q, [TI|, such that at the 
minimum in a, one expects to a good approximation: 

mm a Tl((r) ~ M H . (4) 

Moreover, one can a posteriori check that, at the stability point, where we have an equilibrium 
between the continuum and the non-perturbative contributions, which are both small, the OPE is 
still convergent and validates the SVZ-expansion. The previous approximation can be improved 
by working with the double ratio of moments f] : 

n BB ,(x) = ^*mt, (5) 

provided that each ratio of moments stabilizes at about the same value of a. In this case, there is a 
cancellation of the different leading terms like the heavy quark mass (and its ambiguous definition 
used in some previous literatures), the negligible continuum effect (which is already small in the 
ratio of moments), and each leading QCD corrections. 



2 Test of the 1/m-expansion 

For this purpose, we use the complete horrible results expressed in terms of the pole mass to order 
a s given by |TT| and checked by various authors |J, which we expand in series of 1/m with the 
help of the Mathematica program. By comparing the complete and truncated series in 1/m, one 
can notice that, at the c and b mass scales, the convergence of the 1/m-expansion is quite bad due 
to the increases of the numerical coefficients with the power of 1/m and to the alternate signs of 
the 1/m series. This feature invalidates the analysis in Ref. 



3 Balmer-mass formula 



The Balmer formula derived from a non-relativistic approach (m — > oo) of the Schrodinger levels 
reads for the Sf vector meson (16| (see also |r?l, [pfl ): 



M s s 



static 



2m 



1- -a 2 + 0.23^- — 



9 



ma. 



-(a s G 2 ) 



(6) 



It is instructive to compare this result with the mass formula obtained from the ratio of moments 
within the 1/m-expansion. Using the different QCD corrections, one obtains the mass formula at 
the minimum in o of 1Z: 



M r ~ 2m b ( 1 + — a 2 s ) M. 



static 
Sf , 



at the value: 



V °~ coul 



9 



4ma s y / 7r 



0.85 GeV" 1 . 



(8) 



2 This method has also been used in for studying the mass splittings of the heavy-light quark systems. 



One can recover by identification in the static limit (mj, - > oo) the previous Balmer formula, but 
a new extra a 2 correction due to the v 2 (finite mass) terms in the free part appears here (for some 
derivations of the relativistic correction in the potential approach see [0]), and tends to 
reduce the coulombic interactions. On the other hand, at the 6-quark mass scale, the dominance 
of the gluon condensate contribution indicates that the 6-quark is not enough heavy for this system 
to be coulombic rendering the non-relativistic potential approach to be a crude approximation at 
this scale. The extension of this comparison to the cases of the Sf-S^ hyperfine and S-P splittings 
is not very conclusive, as in the former, one needs an evaluation of the correlator at the next-next- 
to-leading order for a better control of the a 2 terms, while, in the latter there is a discrepancy 
for the coefficients of the gluon condensate in the two approaches, which may reflect the difficulty 
of Bell-Bertlmann JlO|, |Tl| to find a bridge between the field theory a la SVZ (flavour-dependent 
confining potential) and the potential models (flavour-independence). 



4 Leptonic width and wave function 

Using the sum rule Lh arid saturating it by the vector Sf state, we obtain, to a good approximation, 
the sum rule: 

e 25mM v a ,--3/2 

M v T v ^ e+e - ~ (ae Q ) 2 

(9) 



8 . 47T 2 c/2 

1 + -^7rama s — — [a s G )ma 1 
3 9 



where cq is the quark charge in units of e; 5m = My — 2m is the meson-quark mass gap. In the 



case of the 6-quark, we use |T3[ 5m ~ 0.26 GeV, and the value of o~ m in given in Eq. (8). Then: 

r T(5 3)^ e+e - ~ 1.2 keV , (10) 

in agreement with the data 1.3 keV. However, one should remark from Eq. (9), that the a s 
correction is huge and needs an evaluation of the higher order terms (the gluon condensate effect 
is negligible), while the exponential factor effect is large, such that one can reciprocally use the 
data on the width to fix either a s or/and the quark mass []. In the non-relativistic approach, one 
can express the quarkonia leptonic width in terms of its wave function \P(0)q, from which, one 
can deduce: 

16tt|*(0)|J (l - 4C F ^ ~ 18.3 GeV 3 . (11) 

Using the expression of cr cou i, one can find that, to leading order, the two approaches give a similar 
behaviour for ^(0)q in a s and in m and about the same value of this quantity. However, one 
should notice that in the present approach, the QCD coupling a s is evaluated at the scale a as 



dictated by the renormalization group equation obeyed by the Laplace sum rule |L2| but not at 
the resonance mass! 



Gluon condensate from M^ s ^ — M^^u 



The value of er, at which, the S'-wave charmonium ratio of sum rules stabilize is: |plj| : a ~ 
(0.9 ± 0.1) GeV -2 . We use the range of the charm quark pole mass value m c ~ 1.2-1.5 GeV to 



3 Larger value of the heavy quark mass at the two-loop level corresponding to a negative value of 8m, would 
imply a much smaller value of the leptonic width in disagreement with the data. 



order a s accuracy JT3[ 0, and the double ratio of the vector V(Sf) and the pseudoscalar P(Sq) 
moments, which controls the ratio of the meson mass squared. The exact expressions of the 
relativistic, Coulombic and gluon condensate lead respectively to corrections about 0.5, 2 and 7 
% of the leading order one. One can understand from the approximate expressions in 1/m that 
the leading corrections appearing in the ratio of moments cancel in the double ratio, while the 
remaining ones are relatively small. However, the expansion is not convergent for the a s -term 
at the charm mass. Using for 4 flavours a g (cr) ~ 0.48tg;}o! an d the experimental data [23]: 
IZyp = 1.082, one can deduce the value of the gluon condensate: 

(a s G 2 ) ~ (0.10 ± 0.04) GeV 4 . (12) 

We have estimated the error due to higher order effects by replacing the coefficient of a s with the 
one obtained from the effective Coulombic potential, which tends to reduce the estimate to 0.07 
GeV 4 . We have tested the convergence of the QCD series in a, by using the numerical estimate 
of the dimension-six gluon condensate g(f a b c G a G b G c ) contributions given in [20). This effect is 
about 0.1% of the zeroth order term and does not influence the previous estimate in Eq. (12), 
which also indicates the good convergence of the ratio of exponential moments already at the 
charm mass scale in contrast with the q 2 = moments studied in Ref. [fTl We also expect 
that in the double ratio of moments used here, the radiative corrections to the gluon condensate 
effects (available in the literature |2f|) are much smaller than in the individual moments , such 
that they will give a negligible effect in the estimate of the gluon condensate. 



6 Charmonium P-wave splittings 



The analysis of the different ratios of moments for the P-wave charmonium shows [PH , |Tl~] , |20"|| that 



they are optimized for: a ~ (0.6 ± 0.1) GeV 2 . The predictions for the scalar Pq - axial Pf and 
the tensor P|-axial Pf mass splittings, given in SN, are satisfactory within our approximation. 



7 a s from the P\ - Pf axial mass splitting 

The corresponding double ratio of moments has the nice feature to be independent of the gluon 
condensate to leading order in a s and reads: 



ML 



M 



~ 1 + a s A* 3 (exact) = 0.01 1 



*+0.008 



(13) 



The recent experimental value 3526.1 MeV of the P\ mass denoted by h c (lP) in the PDG compi- 
lation j24j almost coincides with the one of the center of mass energy, as expected from the short 
range nature of the spin-spin force Q. Using a naive exponential resummation of the higher order 
a s terms, we deduce: 

a^a' 1 ~ 1.3 GeV) ~ 0.64lg;f| ± 0.02 (14) 

which implies: 

a s (M z ) ~ 0.127 ±0.009 ±0.002, (15) 



4 For a recent review on the heavy quark masses, see e.g. 23 
5 Seee.g. M. 



where the error is much bigger than the one from LEP and r decay data, but its value is perfectly 
consistent with the latter. The theoretical error is mainly due to the uncertainty in A a , while 
a naive resummation of the higher order a s terms leads the second error. Though inaccurate, 
this value of a s is interesting for an alternative derivation of this fundamental quantity at lower 
energies, as it can serve for testing its g 2 -evolution until Mz- Reciprocally, using the value of a s 
from LEP and r-decay data as input, one can deduce the prediction of the center of mass (c.o.m) 
of the P] states given in Table 2 of SN. 



8 T — r] mass splitting 

For the bottomium, the analysis of the ratios of moments for the S and P waves shows that they 
are optimized at the same value of a, namely [PH| : a = (0.35 ± 0.05) GeV~ 2 , which implies for 
5 flavours: a s (a) ~ 0.32 ± 0.06. We shall use the conservative values of the two-loop 6-quark 
pole mass: m& ~ 4.2 — 4.7 GeV. The splitting between the vector T(Sf) and the pseudoscalar 
77&(Sq) can be done in a similar way than the charmonium one. One should also notice that, to 
this approximation, the gluon condensate gives still the dominant effect at the 6-mass scale (0.2% 
of the leading order) compared to the one .08% from the a s -term. However, the 1/m series of 
the QCD a s correction is badly convergent, showing that the static limit approximation can be 
quantitavely inaccurate in this channel. Therefore, one expects that the corresponding prediction 
of (13+Io) M eV is a very crude estimate. In order to control the effect of the unknown higher 
order terms, it is legitimate to introduce into the sum rule, the coefficient of the Coulombic effect 
from the QCD potential [pH], which leads to the "improved" final estimate: 



Mr - M Vb « (63-H) MeV, (16) 
implying the possible observation of the rjb from the T radiative decay. 



9 Gluon condensate from — M Xb 

As the S and P wave ratios of moments are optimized at the same value of a, we can compare 
directly, with a good accuracy, the different P states with the T (Sf) one. As the coefficients 
of the a 2 corrections, after inserting the expression of o~ C oui, are comparable with the one from 
the Coulombic potential, we expect that the prediction of this splitting is more accurate than in 
the case of the hyperfine. The different double ratios of moments leads to the predictions of the 
Xb states Pq , Pf , P 2 3 given in Table 2 of SN in good agreement with the data, if one uses the 



value of a s |19| and of the gluon condensate obtained previously. Reciprocally, one can use the 
data for re-extracting independently the value of the gluon condensate. As usually observed in the 
literature, the prediction is more accurate for the center of mass energy, than for the individual 
mass. The corresponding numerical sum rule for: 

< M c.o.m _ M ^ (17) 

is due to + (1.53^42) x 10~ 2 of the relativstic effects, + (l.20loi) x of the Coulombic and 
+ (0.28±[^)GeV- 4 (a s G 2 ) of the gluon condensate ones. It leads to: 

(a s G 2 ) ~ (6.9 ± 2.5) x 10~ 2 GeV 4 . (18) 

We expect that this result is more reliable than the one obtained from the — M Vc as the latter 
can be more affected by the non calculated next-next-to-leading perturbative radiative corrections 
than the former. 



10 Update average value of (a s G 2 ) 



Considering the most recent estimate (7± 1) x 10~ 2 GeV 4 from e + e~ — > J = 1 hadrons data using 
r-like decay [[| as an update of the different estimates from the light quark systems (see Table 2 
of SN), we can deduce from Eqs. (12) and (18), the update average from a global fit of the light 
and heavy quark systems: 

(a s G 2 ) ~ (7.1 ± 0.9) x 1(T 2 GeV 4 . (19) 

This result confirms the claim of Bell-Bertlmann Jl(J |TTJ stating that the SVZ value |l[ has 
been underestimated by about a factor 2 (see also |57|, |28[).More accurate measurements of this 
quantity than the already available results from r-decay data || are needed for testing the previous 
estimates from the sum rules. 



11 Conclusions 

We have used new double ratios of exponential sum rules for directly extracting the mass-splittings 
of different heavy quarkonia states, the value of the gluon condensate and of the QCD coupling 
a s . Our numerical results are summarized in Eqs. (12), (18) and (19) and in Table 2 of SN, 
where in the latter a comparison with different estimates and experimental data is done. We 
have also succeeded to derive the non-relativistic Balmer formula from the sum rule, using al/m- 
expansion, where we have also included new relativistic corrections due to finite value of the quark 
mass. However, this expansion does not converge at the c and b quark mass scale. 
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